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INTRODUCTION 

The development of heavy transport and military aircraft with 
high performance characteristics has placed increased demand on 
landing gear and tire components. Specifically, aircraft tires 
have developed through the years to a point where the ratio of tire 
weight to aircraft weight has steadily decreased, due in part to 
improved tire cord materials, rubber compounds, and tire struc- 
tural design. However, aircraft tire operations are more demand- 
ing wherein the ultimate strength limits of the tire carcass are 
being approached. Thus, a need exists for the development of 
methods for predicting the strength limit of aircraft tires on 
some rational basis other than indoor dynamometer testing. 

The problem of tire strength differs from that encountered 
in the strength of normal metallic structures and the reasons lie 
in the complexity of the material characteristics of the constitu- 
ents in an aircraft tire. In an aircraft tire, almost all loads 
are carried directly by the textile cord reinforcement, commonly 
nylon as is presently practiced. Some compressive loads are car- 
ried by rubber elements in the tread but these are rather minor. 

In addition to the strength problem, aircraft tires generate sub- 
stantial temperatures as they roll, both during takeoff and 
landing. The temperatures are not high in the usual sense of 
metallic materials, being on the order of hundreds of degrees 
Farenheit or less, but they are high enough to seriously degrade 
the strength of the load-carrying textile components in the tire. 
Excessive temperatures can easily cause a tire to fail even if the 
stress levels are not overly high. 

In addition to the interaction of temperature and stress 
even more complex phenomena are at work in the repeated daily 
use of aircraft tires. These phenomena are basically associated 
with the long term strength and adhesion degradation of an air- 
craft tire due to cyclic variations in stress and temperature. 

This degradation is an area where not much information is avail- 
able, and represents a region where the field usage of tires must 
depend on experience and judgement rather than upon science. 
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Because of this strength degradation factor, the useful life of 
aircraft tires is limited, but these limits are not clearly defined 
and they vary from one aircraft to another and from one user to 
another. 

Substantial advances are being made in the tire industry in 
calculating stresses imposed on textile cord structures due to 
load-carrying and inflation pressure. However the area of tire 
temperature build-up, which is a major cause of tire failure, has 
been almost untouched and there is no analytical technique suit- 
able for predicting aircraft tire internal temperatures on a 
rational basis. Hence, the tire designer or manufacturer is not 
able to compare the anticipated temperature build-up in a new 
design with previous designs or with known temperature- strength 
characteristics of the materials which he wishes to use. This 
void is a serious short-coming in the design process, particularly 
in the preliminary design of high performance aircraft tires, 
where a knowledge of anticipated higher temperatures during normal 
take off and landing cycles would be valuable in assessing the 
validity of proposed tire, aircraft weight, and ground speed and 
taxi length combinations. 

The present research effort is directed toward construction 
of a rational method for evaluating internal temperatures of an 
aircraft tire and for calculating these rapidly and easily, so 
that the method may be used as a design tool in both the tire 
and airframe industries. 

GEOMETRY OF TIRE DEFORMATION 

The basic material property used in the development of the 
tire heating computational method is the hysteretic loss charac- 
teristic of polymeric materials during cyclic stressing. This 
characteristic is illustrated in figure 1, where a typical stress 
stjrain curve for a lossy material is illustrated. As the material 
is stressed from point A to point B and return, the figure shows 
that the stress-strain curve of the loading cycle does not coin- 
cide with the stress— strain curve of the unloading cycle. The 
net area under the curve which remains as the material is return- 
ed to its original stress state represents a loss of 
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Fig. 1. Typical stress-strain curve of a material 
illustrating hysteretic loss 


mechanical work during the cycle from A to B and return. This 
lost work is converted into heat and is either diffused through 
the element and eventually lost from its surfaces or must be 
accounted for in some way in a heat balance. From the analytical 
point of view, the work done in a stress cycle as experienced 
by a material point in passing through the contact patch is 
expressed most usefully by the use of a simple linear visco- 
elastic model such as the Kelvin-Voigt solid, and using such 
a model the energy loss per load cycle may be written as in 
Equation 1. 

AE = U TT AV • tan6 (1) 

where 

AE = Energy loss/cycle 
U = Change in Strain Energy from A to B 
tan6 = Loss Tangent of Material 
AV = Volume Element 

The question of whether the energy loss depends on the 
change in strain energy between two non-zero stress points, or 
upon the difference in strain energy levels each calculated from 
zero stress level, is treated in Appendix A. it can be 
demonstrated that the former formulation is the correct one. 
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This formulation presumes that the initial stress induced 
by tire inflation has no influence in a linear sense on the 
energy loss per unit stress cycle. However, the initial stress 
state may influence the material characteristics, which are not 
linear in this case. Methods for calculating the change in 
strain energy will be discussed in the next section of this 
report. 

The method used for calculating the change in 
strain energy is based on two assumptions. The first is that 
the tire is made up of a series of material points, each of 
which undergoes a change in stress state corresponding to the 
points A and B in figure 1, as it moves from point A to B shown 
as in figure 2. 



Fig. 2, Extreme points of stress excursion in a tire 

As a material point moves from A to B in figure 2, it is 
assumed that the general wave form of the stress cycle depected 
in figure 1 is unimportant, and high harmonics of the cyclic stress 
are neglected in favor of a single fundamental cyclic stress change, 
which will be described in terms of an assLimed geometry of the tire 
at point A and B. In other words, no allowance is made in this 
simplified analysis for a Fourier decomposition of the stress 
state as the tire moves from A to B, but rather the fundamental 
harmonic is considered as defining the change in strain energy 
between the two points. 
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The second major assumption made in calculating the strain 
energy of the tire as it rolls through the contact patch is that 
at both points A and B of fig. 2 , the cross-section of the tire 
may be represented geometrically by its neutral axis, around 
which act both membrane and bending strains. 

The neutral axis of the tire section is usually taken at 
the mid- thickness location which may be determined approximately 
by estimation or alternately by calculation. For the case 
of the inflated tire cross section, the neutral axis may be 
located by; 

(a) Use of a tire section to approximate the neutral axis, 
assigning numerical values of co-ordinates to selected points. 

(b) Use of the uninflateid tire dimensions and cord angles 
to provide input data to calculate the inflated tire shape, 
using such methods as the Hof ferberth equation [4 ] , or 
finite element methods. 

In the case of the tires studied in this report, both methods 
were examined. While the Hof ferberth equations are convenient 
and seem to agree well with measurement in the crown and upper 
sidewall, they do not agree well in the lower sidewall and bead 
area where greater bead and rim constraint is present. For that 
reason the numerical data generated for this report were obtained 
by approximating the inflated tire cross-section from drawings. 

For the deflected tire a similar but more severe situation 
exists because the deflected cross section must be obtained. 
Fortunately most aircraft tires are molded in the tread and 
shoulder regions so that a clear line of demarcation is present 
between portions in contact with the runway and portions not 
in contact. This demarcation allows one to assume that the part 
of the cross section in contact is flat when at the center of the 
contact patch, as shown in fig. 3. The sidewall portion, out 
of contact, can then be defined geometrically by a separate 
calculation. Three possible methods for such a calculation include 

(a) Some form of finite element solution involving large 
deformations 

(b) A linearized, energy based deflection analysis 

of the sidewall arch. The boundary conditions imposed 
on the sidewall section are; 
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1. At the shoulder the sidewall must have the same 
deflection as the tread region. 

2, At the bead the arch acts as a built-in section 
with no change in slope. 

(c) The assumption that the high pressures commonly used in 
aircraft tires cause the sidewall to act essentially as a 
membrane, neglecting bending stiffness. For this assumption, 
the sidewall may be approximated by using a segment of a 
circular arc as the sidewall shape. 

For simplicity and computational ease the latter method 
has been used to obtain sidewall shapes. It is surprisingly 
realistic in spite of its linear approximate character. Tire 
cross section widths calculated by this method agree well 
with those observed experimentally. 

Figure 3 illustrates the neutral axis of the structural 
part of the carcass by the lines labeled 0-1-2. The inflated 
tire geometry is presumed to be represented by fig. 3a. The 
deformed tire, under deflection A as shown in fig. 3b, is 
presumed to operate under conditions of being perfectly flat 


0 



(a) 

Assumed geometry at 
point A, Fig 2 

Fig . 3 


1 


1 0 ^ 



Assumed geometry at 
point B, Fig. 2 
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in contact with the runway in section 0—1/ and having a reduced 
radii in the sidewall as given by the symbol r^. The bead, 
point 2, is fixed. The regions 0-1 and 1-2 are considered ihex- 
tensible for computational convenience which allows the complete 
geometric description of the various radii of curvature, onee the 
tire deflection A and the initial tire geometry are given under 
inflated conditions. The computational details of the geometric 
representation are given in Appendix B. 

A material point is subject to both bending and membrane 
strain during deformation of the tire from point A to point B 
as shown in fig. 2. The Ig^ending strains are dependent on the 
radius of curvature of the tire carcass midline as illustrated 
in fig. 3, as well as upon the distance measured perpendicular 
from this carcass midline. The bending strain can be calculated 
in the usual linear elastic fashion based upon a knowledge of 
change in radius of curvature of the carcass midline between the 
two points in question in fig. 2. 

Membrane strains also exist in the tire cross section due 
to the high inflation pressure commonly used in aircraft tires. 
Some of these membrane strains are not directly obtainable from 
the geometric changes which have been postulated, since they 
depend heavily on contact patch boundary conditions. Other mem- 
brane strains must be obtained from a knowledge of the membrane 
stresses in the tire. These membrane stresses in turn depend 
upon the curvature characteristics of the tire, as well as upon 
certain geometric variables which may be calculated from the tire 
cross section geometries given in fig. 3. The details of these 
calculations are given in Appendix B. In particular, the sidewall 
area must be carefully defined since the point of vertical tangen- 
cy of the tire is an important factor in computation of membrane 
stresses. This point of vertical tangency, measured as the 
distance b in fig. 3, shifts during the deformation process, 
and the position of the various material points relative to 
this point of vertical tangency must be tracked in the subsequent 
computations . 
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Since aircraft tire deflections are normally large, a fur- 
ther assumption is made that the tire deformation process causes 
changes in the tire cross section, such as shown in fig. 3, but 
that deformations in the circumferential direction are negli-^ 
gible . 

It is also necessary to track the radius of a material 
point from the axis of rotation of the wheel during the deforma- 
tion process, as shown in fig, 4. The variable used for such 
measurement is denoted by the symbol z» In fig. 4, the point j 
represents a material point which changes its axisymmetric radius 
during the deformation process as the wheel rotates from condition A 
to condition B. The details of the computa.tion of such a change 
are also given in Appendix B. 

Bead tensions also change during the deformation process 
and must be incorporated into a cyclic stress mechanism 
involving a loss characteristic of the bead itself. These bead 



Fig, 4. Tire cross sections 

tensions depend on angle of tangency of the side wall at the 
bead, denoted as point 2 in fig. 3. Computations of bead 
tension are given in Appendix G where the bead geometry 
and assumptions governing it are discussed in more detail. 
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TIRE DEFORMATION AND STRESS 

For the calculation of the strain energy change during 
the rotation of a material point through the contact patch we 
rely heavily on the assumed tire deformation given in figures 
2 through 4 , 



Fig, 5, Tire coordinate notation 

Fig. 5 illustrates the geometry involved and the notation being 
used here. z represents radial distance from the axis of 
rotation of the tire, <p represents the meridional direction and 
0 represents circumferential direction. 

For purposes of computing changes in strain energy states 
during the cyclic loading process the tire is divided into five 
regions, denoted as I-V in fig. 6. 

The rationale behind the division of the tire into these 
parts is that region II represents the tread region of the tire, 
which is considered to be rubber-like and in contact with the 
runway surface. It is usually characterized by material isotropy 
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Edge of 



Fig. 6. Tire cross section showing regions of 
diffferent characteristic and loading 
patterns 

and by an elastic shear modulus which is nearly linear up to 
reasonably large strains. Further, it exhibits near incom- 
pressibility, and that fact is used in computing the strain 
state associated with contact against the runway surface. 

The region denoted by I is immediately beneath the tread 
region of the tire and is firmly affixed to it, so that the 
membrane strain state of the two elements is considered to 
be the same. However, region I is assumed to be plane ortho- 
tropic in its characteristics, and may be thought of as a 
sheet of material lying in the (j>-0 plane as illustrated in 
f ig . 7 . 
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Fig. 7. Laminate characteristics in Region I 

In fig. 7 the diagonal lines represent textile reinforcing cords, 
which serve to give the primary stiffness to the carcass region. 
Because of their diagonal character, which is typical of a bias 
aircraft tire, the element shown in fig. 7 will exhibit different 
elastic constants in the cj> and 0 directions. Such materials have 
been studied in the tire industry and their characteristics can 
be estimated. Reference [6] gives a detailed derivation of the 
appropriate elastic constants for a plane orthotropic material 
using the structure shown in fig 7, adequate to within a range 
of diagonal reinforcing angles a encountered in aircraft tire 
construction. Fig. 7 can be characterized in terms of elastic 
constants dependent upon the rubber and cord reinforcing materials 
and upon the half-angle a illustrated in the figure. 

Regions III and IV of the tire, as illustrated in fig 6, 
represent the sidewall areas which also are subjected to changes 
of curvature as the tire deforms. These areas do not contact 
the runway surface, so that in Region III the carcass strains 
are again determined and the elastic constants obtained from the 
model of fig. 7. These constants are utilized to determine the 
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appropriate strain energy expressions. Region IV represents 
a rubber covering used primarily as a scuff or abrasion resis- 
tant coating on the outside of the tire. It contributes to tire 
mechanical loss and to the subsequent heat build up, but is not 
structural in nature. The assumption is made that this covering 
is bonded tightly to the carcass of Region III and undergoes the 
same membrane strains as those of Region III, but, in addition, 
undergoes bending strains due to its distance away from the 
neutral axis, which because of the relatively high textile modulus 
is usually located in Region III of the carcass. 

Region V is the tire bead. It requires a separate analysis, 
not only because of the force fluctuations in the beads but also 
because the beads are conventionally made of rubber coated steel 
wire which has significantly different elastic and loss charac- 
teristics than the remainder of the tire. 

For computational purposes the calculation of strain energy 
stored in the tire during the deformation process is accomplished 
by means of splitting each region of the tire into a number of 
cells, each cell being denoted by a pair of numbers in matrix 
notation fashion. The element numbering system for the cells in 
the tire cross-section is illustrated in figure 8. The elements 
shown in this figure are purely for illustration purposes. Any 
number of elements may be chosen and the order of them may be 
adjusted to suit the needs of the user, but for the instances to 
be illustrated here numbering begins in the bead area and termi- 
nates at the center line of the tire crown. Numbering of the 
vertical location of the elements begins at the inner surface 
and terminates at the tread or outer surface of the tire. 

The inner liner of the tire, while pure rubber, is thin 
enough so that it may be included with the composite carcass 
material . 

Computation with the individual elements is facilitated by 
use of the matrix notation illustrated in figure 8. However, for 
a general theoretical treatment of the strain energy change in 
the tire in this report the notation of figure 6 is adopted, where 
the major portions of the tire are delineated. An approximate 
analysis of the strain energy in each of the Regions I-V of the 
tire is now given. 
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Fig. S. Element numbering system 


Region I 


In this section the assumed geometry of figure 3 is used to 
develop membrane stresses in the (j) (meridional) direction. A 
shear effect is included and will be discussed later. 

(a) Undeformed Position 

First consider the tire in the undeflected position (point 
A of fig. 3). Referring to figs. 3, 4, and 9, equilibrium 
of the axisymmetric annular shell element o-s gives a 
membrane stress 

P (z^-b^) 

= —2. s 

^<j) 2h z sin (j> 

m s ^s 


Assuming that changes in shape from the uninflated to 
the inflated state are small, the bending stresses here 
are negligible so that 


b 


(3) 
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9. Membrane force nomenclature 


The membrane stresses in the 0 direction are obtained from 
the nature of the load carrying cord network shown in fig. 7. 

It may be demonstrated that in a soft matrix where the cords 
carry the entire load. 


a 


0 


m 



2 

cot a 


(4) 


where a is the cord half-angle shown in fig. 7. The cord angle 
relation is discussed in more detail in Appendix D, The 6 
directed bending stresses a,re also assumed to vanish so that, 
similar to equation (3) 


The total stress state in the undeformed but inflated 
condition, point A of fig. 10, is given by the membrane stresses 
just determined. It is denoted by the subscript 'A' , taken 
from Fig. 10. 
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(a 




Ref. equation (2) 


( 0 ^ 9 )^ ^ ^0,m equation (4) 

(b) Deformed Position 

Eq. (2) gives the meridional membrane stress o_, in terms 

(Pm 

of the sidewall geometry of the tire. In the previous 

sections it was shown how sidewall geometry may be 

determined in the deformed position. Assuming that it is 

available, then oa may be found at the points labeled 

C and E in fig. 11, i.e. the edges of the contact patch. 

Denote this by a , 

(f),m 




Fig. 10 


Fig. 10. 


Side view of tire 
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The entire contact patch region C-D-E-F of figs. 10 and ll 
is assumed to be flat and in contact with the runway. Were it 
not for friction between tire and runway, the flat contact patch 
region would be subjected to a constant membrane stress value 
across its width of 



where a, is the value of the membrane stress at the edge of 
(j) ,m 

contact, points C and E of fig. 11. 

Due to the presence of frictional effects the stress a. _ 
cannot decay fully to the value given by equation (6a) , but instead 
a correction factor must be placed on the stress difference 
depending on position inside the contact patch. This analysis 
is given in Appendix K, and the appropriate factors are included 
in the numerical computations which follow. 

The circumferential membrane stress may once again be 
determined from the relation originally given in equation 4. 


a 


0 , m 



cot a 


(6b) 


With this relation, the membrane stress state may be determined in 
the deformed position of Region I. This stress must also be 
corrected by the same friction effect. 

In addition to these membrane stresses, bending effects 
are present in the contact region of the tire since the carcass 
is flat when it is in contact with the runway surface. The 
curvature changes undergone by the carcass are assumed to be 
exactly equal to the reciprocals of the principal radii of 
curvature of the tire, namely the values R of Fig. 10 and r^ of 
Fig. 3. Using these values, the bending strains may be determined 
as 


S,b = = (?-y)(i+§) 

^ ' c 

where y is measured from the inside surface positive toward the 
outer surface of the tire. Note that when running on a drum 
this curvature can be included by defining the curvature change 
in the e direction as (1/R+2/D) , where D is the drum diameter. 
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These strains may be used in conjunction with the usual 
orthotropic stress-strain relations commonly applied to bias 
cord-rubber constructions to obtain the stresses associated with 
bending. These become 


/b = 






E, 


0,b 1 u u '•^0.b'^^0<j)^(|.,b^ 

(J)0 0c|) 


(7) 


(8) 


The total stresses are now given by the sum of the bending 
terms equations (7) and (8) added to the membrane stress terms 
of equation (6). Again^ the subscript 'B' is from fig. 10 to 
denote the deformed condition 




= 

+ a. , 

0,m e,b 


The strain energy density is now based on the stress 
differences which take place between points A and B of fig. 10 
and which were illustrated earlier in fig. 1. 

In order to form this strain energy density the stress 
difference or change from A to B is first formed. This stress 
change is then used to determine the strain change accompanying 
it. These changes then can be combined to form the change in 
strain energy from the A to B conditions. Denoting changes by 
the symbol A, one finds 


Aa, = (a ) - (a ) 

♦ ^ B ‘•'a 


= <%>3 - 


Using orthotropic elastic properties, 

AiJ , A a „ 


a€j 


4 ) 


“ M 


0 E 


( 10 ) 


Aa 

^ 04 ) ~E 
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( 11 ) 




Aa 


£ 

e 


i 


The strain energy density is now given by 


( 12 ) 


The elastic constants needed for these formulations of 
strain energy are derived in Appendix E. 

The other contribution to strain energy in Region I is 
caused by the change in the membrane stress state as the tire 
rolls through the contact patch region. These changing tensile 
stresses result in shear in the sidewall and contact area. This 
analysis is rather lengthy but leads to contributions which can 
be significant, particularly for larger tires. The details 
of this analysis are also presented in Appendix E. The resulting 
strain energy density is given in summary form 


■" f ^ . 9 

2 

0.858{(a, ) ”(cT, ) } cot ads 


-|.,m ^ ^,m 3 Jq 

1 

{3A (D-A) } 



where the symbols used are defined in Appendices E and F. 

Because of the complications arising from introduction of 
anisotropic elastic constants, these two strain energy densities 
will be summed algebraically in order to give the total strain 
energy 

U = U + U (14) 

c s ' 


for Region I. 


Region II 

The basic assumption used in calculating the strain energy 
in the tread rubber is that this material is non-structural , and 
that it undergoes the same membrane strain as the carcass 
immediately underneath to which it is bonded. Its bending strain 
is determined by its distance from the neutral axis and by the 
Kirchoff hypothesis. 

(a) Undeformed Position 

The carcass membrane stresses in the undeformed position have 
been determined in equations (2) and (4). These stresses can be 
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used along with the conventional orthotropic stress- 
strain relations to obtain the membrane strains in that 
position 


(j) ,m 


e“ 


(""e) 


- y 


e E, 


<€e>, = 

A ' A 


(^e) 

eT 


- y 


0(f) eT" 


( 15 ) 


(b) Deformed Position 

Similarly, the membrane stresses in the 
found in Eqs. (6a) and (6b), from which 
strains are found to be 
* * 




(|),m 


^0 


■0 ,m 


0 ,m 

E. 


- y 




0<j) E 


<t> 


deformed carcass were 
the corresponding 
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In addition the deformed state also exhibits bending strains 
Due to different distances from the neutral axis, these are 
now dependent on local element geometry. in general the bending 
strains may be expressed as 


= (y-y)-^ 

C 

( 17 ) 

^0,b (y-y) (r + 5) 


by 


The total strains are now given for the deformed position 


J 

^ B 



(6. ) 

^ B 


^^0,b 


( 18 ) 


The changes in strain between the undeformed 'A* position 
and the deformed 'B' position are taken to be principal strains. 
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a£, = (€,) - (€^) 

♦ T g <P 


a€„ = (6„) - (€J 

® ® B ® A 


(19) 


These strains are caused by deformation of the tire carcass, 
and result in stress changes given by 


Aa 


i-y‘ 


(A 6, + yA£.) 


Aa, 


E 


1-y 


j(a£q + pa€ ) 


where and y refer to isotropic elastic constants of the tread 
rubber. Thus, strain energy can be expressed in the form 


U 


T 


2G[a€^ + A£^ + 


( 20 ) 


The maximum value of the compressive contact stress results 
in compressive strains in the approximate value of 



where the rationale for the use of the numerical value 1.875 
is given in Appendix G, where p^ is the tire inflation pressure 
and E the Young ' s modulus of the tread rubber . 

In spite of the incompressible nature of the tread rubber, 
all three principal strains AEq, and can exist since the 

groove pattern molded into the tread allows bulging of the indi- 
vidual tread ribs independent of the (j) and 0 strains. 

The additional strain energy caused by the contact pressure 
is given by 



- 21 - 


U 

z 


1 

2 


(1.875 p^) ^ 


E 

r 


(21) 


The total strain energy change is now given by the sum of 
equations ( 20 ) and ( 21 ) 


U = 


Um + U 

T z 


2Gl&el + A£^A€g] + I 


(1.875 P^) 
E 


Region III 

(a) Undeformed Position 

In this region the analysis of strain energy density is very 
similar to that of Region I, In the inflated state, as 
shown at point A of fig. 10, the membrane stress in the 
meridional direction <J> is given once more by equation (2) 


P / 2 ,2, 

(n ) - a = — ^ 

(J) (j),m 2 h z sin tj) 

and the circumferential stress by 


equation (4) 




A 


m 


4-2 

a, cot a 
(J> ,m 


Again, bending stresses in both t() and 0 directions are 
considered negligible. 


(b) Deformed Position 

In the deformed position the membrane stresses are still 
governed by equilibrium considerations such as given by 
equations (2) and (4), but now with use of the deformed geometry 
appropriate to the center of the contact patch such as 
point B of fig. 10. This computation still relies on the 
assumed axisymmetric character of the stress, so that 
the tire acts as if it were inflated in a barrel. The 
new sidewall positions must be calculated based on the 
deformed sidewall geometry as discussed earlier. Using 
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this assimied geometry, then the deformed position stresses 
become 


a 


^,m 


il 

2h - 
z 


2 t-2 

~b 

sin 


)_ 

J 


a 


e ,m 


^2 

a, cot a 
(p ,m 


(22) 


where the symbols z, b, (J) denote those quantities measured 
in the deformed position. 

Bending strains are also present in the (f) direction due to 
the large changes in radius of curvature which are imposed on 
the tire sidewall. Curvature changes in the 0 direction are 
neglected leading to a strain change 

b = (y-y) 

f 1 

where r^, = radii of curvature in the meridian plane, final 

and initial respectively. 

Using Hooke's law for an orthotropic material, the bending stresses 
generated by this bending strain are 


(j>,b 


E 


^ ^(j)0^e<j) 




rb 




The total stresses in the deformed state are 


(. 24 ) 


= °e,m + ■’e,b 


(. 25 ) 
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and the stress excursion from position A to position B is 

Aa. = (a ) - (a.) 

« 4' B 4 a 


( 26 ) 


Aa = (a ) - (a ) 

® ® B ® A 


This stress change causes a strain change 

AC 


1 _ II ^ 


A6, 


E 


E. 


4,0 E, 


(27) 


Aa 


- y 


Q<\> E 


These changes can be used to form the strain energy 

+AagA€g) 


(28) 


The shear forces previously discussed are also present in the 
sidewall, and these should also be included here as given in 
Eq. (13 ) , so that the total strain energy is given by 

U = U + U 
c s 


Region IV 

Here the same approach is taken as in the previous Region 
III where the total sidewall strain change is calculated using 
both the membrane effects and the bending induced by the 
curvature changes. These strains are given by equation X28) . 
Using this equation and knowing the isotropic character of the 
rubber covering on the sidewall and shoulder area not in contact 
with the runway surface, then 

U = 2G(A€2+i€2^A£g-Ae^) (29) 

as in equation (20). Shear deformation is neglected here. 
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Region V 

In this region the bead elements themselXres are considered 
as components which lose energy during cyclic stressing of 
the tire. An analysis of this is given in Appendix C, and 
the result of that analysis may be used to introduce a total 
or gross strain energy change in each bead of the tire. Using 
the expressions derived in Appendix C, the energy change becomes 

U = |(Tf-T.)VA^Eg (30) 


where 


final bead tension 
initial bead tension 
bead cross sectional area 
bead modulus 



I 
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HEAT GENERATION RATES 

For purposes of obtaining heat generation rates from strain 
energy stored during elastic cycling of the tire, the conventional 
Kelvin-Voigt mechanical model of an elastic solid is utilized 
which exhibits loss characteristics under cyclic stress. Figure 
1 can represent experimental data observed on many materials, 
including the polymeric materials used in aircraft tires such 
as rubber, nylon in textile cord form, and even the bead wire 
assembly, provided that the loss characteristics are measured 
properly. One widely accepted method of representing such a 
phenomena, the Kelvin-Voigt model, is illustrated in fig. 12 
using the conventional spring and dashpot notation. 


K 



Fig. 12. Spring-dashpot model of a material 
exhibiting loss 

Imposing a cyclic sinusoidal strain in the form 

e = e^sinwt, (31) 

the stress necessary to maintain such a sinusoidal strain is 
given by 

a = Ke + Ce = Ke^sinwt + Cw 0 ^coscot (32) 


I 
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It is customary to convert this representation into a form using 
the elastic constants E' and E", which represent the character- 
istics of a material exhibiting loss, as follows 


a = E*e sinwt + E"e cosot 
o o 


(.33) 


It is clear from this representation that both E' and E" may be, 
and usually are, functions of frequency. In addition, particu- 
larly in polymeric materials, they show marked dependence on 
temperature and, of course, on material composition. 

The energy loss during one cycle of stress is given by 


Energy Loss = W 


ade = ■irE"£^ 

half- 2 

cycle 


Defining the ratio of E" to E' as follows 


(34) 


E"/E' = tans 


(35) 


and combining equations 33 and 34 gives 
E * G 2 

W o tanS — irUtanS 

2 


(36) 


where U is the change in elastic strain energy during the cyclic 
stressing process. Equation 36 is in the form in which energy 
loss is commonly expressed, using the elastic constants E' and 
E" directly obtained from experimental data, or, alternately, 
using tan6 directly obtained from measurement. Several instru- 
ments are commercially available for the measurement of loss 
properties such as tan6, and values of tan6 are available for 
many polymeric materials. 

The energy lost in going from one stress state to another 
is represented by using the change in stress or in strain when 
computing the strain energy term U in equation (36). This energy 
loss has been discussed and is demonstrated experimentally in 
Appendix A. 
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The units of Eq. 36 are the same as those of strain energy, 
namely work per unit volume. However, these units must be converted 
to a rate of heat generation per unit volume in order to be used in 
the general heat diffusion equation discussed in the subsequent 
sections. Since aircraft velocities are commonly expressed in 
English units of feet per second, tan5 is dimensionless, and 
the tire diameter may conveniently be expressed in inches, and 
using the elastic strain energy in inch/pound units, then one 
may write; 


U/12 = — Ibs^ — 6 ne r gyy stress cycle 

' cu. in. volume 


1 3 

X BTU/in /cycle 

3 

X 252 cal/in /cycle 


But 


r CO 

r 

12 


1 3 

X cal/cm /cycle 

(2.54) ^ 

= V 

o 


where v. 


10 


aircraft velocity, fps 
rolling radius of tire, inches 

- A/3 

tire angular velocity, rad/ sec 
12v . 

-I r-T^ X cycles/sec. 

(r^-A/3) 2 tt ^ 


Loss/cycle = W = irUtanS 

. 3 

Heat generation rate q in cal/cm /sec may be expressed as follows 


- U X ^2 ^ 778 ^ ^ 


2.54 


X 


12 

2ir 


(r^-A/3) 


X TTtan6 


= .00988 


Uv tan5 
o 

(r^-A/3) 


Uv. 


0 . 01 ' 


tan(S 


( 37 ) 
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where r = outside radius of tire, inches, 
o 

A = tire deflection, inches 

= aircraft velocity, feet/second 

= elastic energy, inch-pounds/cubic inch, stored 
due to the change in stress state. 


U 
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TEMPERATURE DISTRIBUTION MODEL 

Referring first to fig. 8, the tire must be divided into 
segments or cells wherein in each the rate of heat generation 
is considered constant, and the temperature can be 
calculated on the basis of being averaged across the cell area. 

In assigning the geometric location of cells such as illustrated 
in fig. 8, material constancy is maintained across each 
cell, in the sense that the interface between the pure rubber 
region and the carcass region is used as a dividing line for 
a series of cells. This procedure allows the properties of each 
cell to be essentially uniform, which is a considerable con- 
venience in subsequent calculations. The stress state in each 
cell is calculated on the basis of the centroid of each cell 
area, and is considered to be averaged across the cell volume. 

The geometric location of the cell is taken to be the coordin- 
ates of the cell centroid of area such as shown in fig. 13 . 



Fig. 13, Distribution of computational elements 
in the tread region 
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Quadralateral elements have been consistently utilized in 
constructing the numerical analog to this solution. Thus each 
element either loses or gains heat from four faces, and this 
exchange may be through conduction when the element in question 
is bounded by a similar element or it may be by convection from 
a surface face, or it may be by conduction to either the flange 
of the wheel or to the runway surface for those special elements 
which come in contact with them. 

For purposes of conduction to the wheel flange, a similar 
set of elements to those shown in fig. 13 must be laid out on 
the wheel flange using the appropriate thermal characteristics 
of the wheel material, so that heat flow can take place by 
direct conduction between the bead area of the tire and the 
flange. There will inevitably be a certain contact resistance 
between the bead seat and the flange material, but this factor 
is as yet unknown and is neglected in the subsequent computations. 

Contact with the runway surface is more difficult to 
quantify since it occurs only periodically. To a first approxi- 
mation that process is treated in a similar fashion using the 
thermal conductivity of the runway material, but allowing only 
that fraction of the conduction process to occur as is propor- 
tional to the time spent by each surface element in contact with 
the runway. The remaining part of the cycle is occupied by 
those tread surface elements which contacted the runway surface 
now freely moving in the atmosphere, and during that fraction of 
the cycle the normal convection process from the surface occurs. 

Using the well known expression for diffusion of heat in a 
solid body under conditions of internal heat generation, such as 
given by Carslaw and Jaeger [1], the equation for conduction 
of heat for a typical element is given by 


KV^e = P If - q (38) 

This expression may be written for finite difference purposes 
in the form 


= q 


'** “k-1,K, £ ^®K-1, 


“k,K+1, £ ^®K+1, £ ®K,£^ 


A0 

At 
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"k, £-1, “k, £, £+l^®K, 

( 39 ) 

where the elements used in equation 39 are illustrated in fig. 
14, and where 

- KA 

“m, n d 

with 

K = thermal conductivity 

A = area of contact between subscripted elements 
d = distance between centroids of subscripted elements 




K.X-I 




K-1.X 

KX 

K*[,l 













Fig. 14. Notation for thermal elements 

The basic method of computation is to use small time 
increments and to calculate the temperature rise in each cell 
utilizing the temperatures of adjoining elements from the 
previous time increment computation. There are several factors 
which make this type of computation possible. First of all, 
it is well known in heat transfer calculations that if the 
quantity given by the expression 
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M 


2 


AX 

aAt 


where 

a = q/ci 

Ax = element size 
At = time increment 


(40) 


is large, then stability of such a time based forward computa- 
tion is assured and convergence to correct results will be ob- 
tained. For that reason cell elements have been kept moderately 
large and time increments small since this reduces requirements 
for computer memory and also insures stability of the computation. 
Secondly, the dominant term is the rate of heat generation q, and 
for purposes for aircraft taxi and take-off the heat Duild up 
occurs over a relatively short period of time, so that the 
heat generation term almost completely controls the temperature 
rise in the tire. Heat diffusion, at least for the take-off 
case, is a minimal effect. Thus, heat conduction 

using adjoining cell temperatures from the previous time interval 
incurs only very small errors in the final temperature distri- 
bution in the tire. 

The heat transfer coefficients from tire to air are not 
particularly well known in this case, but again it is fortunate 
that taxi take-off cycles are relatively short, so that again 
the rate of internal heat generation is the dominant term. 

However, in studying taxi and combined taxi take-off operations, 
more emphasis is placed on the accuracy of the heat transfer and 
heat conduction coefficients associated with this problem. In 
that case those coefficients become more important. The numerical 
proportioning of the heat transfer coefficients is given in 
detail in Appendix H. 

The method of temperature computation using Eq. 39 now 
consists of starting from whatever initial temperature state 
one wishes in each of the cells, allowing a short time interval 
to take place and computing the heat generated in each cell using 
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the concepts developed in the previous two sections. In 
carrying out such computations a knowledge of the material loss 
characteristic, tan6, must be known, and these are given in the 
last section of this report. 

At the end of the first time interval the temperature field 
in the entire tire is calculated and is used as a basis for 
the heat transfer characteristics during the second time 
interval, which follows immediately. The computation is 
sequential, thus requiring minimum memory storage in a computer, 
and may conveniently be carried out on relatively unsophisti- 
cated equipment such as minicomputers. 

The contained air temperature in the tire is often measured 
and used as an indication of overall tire heating. Its 
calculation requires a separate analysis given in Appendix I . 
From this analysis, the contained air temperature is shown to 
lag behind the inner surface temperature of the tire. In the 
initial stages of tire heating, its value is given by 

Ay = 3At(e-y) (41) 

where 

y = contained air temperature 

0 = weighted average of tire inner surface temperatures 

At = time increment 

3 = a constant involving properties of the air and heat 

transfer coefficients. 

The heat transfer coefficients chosen for the inner surface 
are presently based on the concepts that inner peristaltic 
pimping velocities are small, just sufficient to insure good 
mixing. This is an area where more measurements are needed. 
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EXPERIMENTAL PROGRAM 

The problem in hand is an extremely complex one and could 
not be expected to be verified without considerable experimental 
data. A number of assumptions are used in the computation of the 
rate of internal heat generation in the tire, and many of these 
assumptions can only be verified by direct measurement of temper- 
ature fields in actual tires. For that reason a comprehensive 
experimental program was carried out by the Impact Dynamics 
Branch, Langley Research Center, National Aeronautics and Space 
Administration, Hampton, Virginia. Most of the work was done 
at the Wallops Flight Cehter, using a specially designed 
truck with a yoke and loading apparatus designed for the testing 
of tires. 

Thermocouples were installed at the University of Michigan 
in approximately twenty 22 x 5.5 8-ply and 12-ply rated aircraft 
tires furnished by NASA Langley. Most of these thermocouple 

installations were done by implanting appropriate thermocouples 
in holes drilled in the tire carcasses after having been buffed 
but prior to retreading. Subsequently the tires were retreaded 
by a commercial retreader, and were wired for internal tempera- 
ture measurements during both straight line rolling as well as 
cornering and braking operations. For purposes of this report 
only the straight line rolling data will be used. Due to load 
limitations of the test equipment the full rated loads of the 
tires could not be applied, and since it was desired to achieve 
normal tire deflections, inflation pressures were also kept below 
normal values. 

Thermocouple locations are illustrated in figures 15 and 16 
for the 22 x 5.5 tires. 

The tire temperature data obtained from typical test runs 
is presented in the next section of this report where it is com- 
pared with computation. In general the use of appropriate 
thermocouples carefully implanted in the tire was sufficient to 
give a good picture of the internal temperature distribution. 

A thorough summary of this information is to be presented in a 
separate NASA document. 

Due to limited load and pressure combinations which could 
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Fig. 15. Element sectioning for thermal analysis. 

Large solid circular symbols represent 
approximate thermocouple - locations in the 
size 22x5. 5| 8 PR test tires 
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be used on the NASA equipment, the University of Michigan con- 
ducted a series of taxi-takeoff roll tests on instrumented tires 
of the same size at Wright-Patterson Air Force Base, Dayton, 

Ohio, using their 120-inch aircraft tire dynamometer. With this 
equipment, full rated loads and speeds could be obtained. These 
data are also presented in the subsequent section and compared 
with computed values. These data represent a different set of 
temperature rise values since they involve loads, pressure^ and 
speeds which are consistent with real operating conditions. 

Referring to earlier sections of the report, it is clear 
that the rate of heat generation depends on the strain energy 
stored in the tire, which in turn depends both on bending and 
membrane effects. While deflection itself controls the bending 
contribution to this problem, the membrane contribution is clearly 
dependent on the inflation and thus both effects are needed to 
give a realistic description of tire temperature rise during ser- 
vice conditions. 

The 22 X 5.5 tire is a relatively small aircraft tire, and 
since one of the most important applications of this type of cal- 
culation is to larger tires, it was felt that experimental 
data should also be obtained on them. Three 40 x 14 aircraft 
tires were also made available by NASA, and these were again 
fully instrumented with thermocouples at the University of Michi- 
gan. Approximately forty thermocouples were installed in each 
tire, using locations near the outer surface, near the carcass 
midline^ and near the inner surface. The points of installation 
ranged from the crown centerline to very close to the bead. 
Arrangements were made to have these tires run on the large dyna- 
mometer facilities at Wright Patterson Air Force Base under the 
Air Force Flight Dynamics Laboratory. Full rated loads and tire 
deflections could be obtained there, but the tests were limited 
in speed in order to obtain longer thermocouple lead wire life. 

A complete series of tests were run on these three tires, in- 
volving various deflections at constant pressure and varying 
pressure at constant deflection. Thermocouple locations for 
these tires are illustrated in figures 17 and 18. 




Fig. 17. Element sectioning for thermal analysis. 

Solid circular symbols represent approxi- 
mate thermocouple locations in the size 
40x14, 22 PR test tires 








Fig, 18, Element sectioning for thermal analysis. 
Solid circular symbols represent approxi- 
mate thermocouple locations in the size 
40x14, 28 PR test tires 
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Tires of two ply ratings, 22 and 28, were used in this 
second series of experiments. Due to considerable differences 
in thickness of the carcass of these two tires ^ they are 
treated in this report as two separate tires. 

While a great deal of temperature data were obtained by these 
two series of test programs, it is so voluminous that only selected 
parts have been used for purposes of comparison with calculations 

in this report. In order to provide a complete record of 
the test data available, however, Table 1 shows the load and speed 
combinations available in straight line rolling from the 
22 X 5.5 NASA tests, while Table 2 -shows the load, pressure 

and speed combinations available from the University of Michigan 
tests carried out on that tire size at Wright-PattersonAir Force Base. 

Table 3 gives a similar summary of test conditions for the 
40 X 14 tires. 

Some of these data will now be extracted in detail and used 
in a subsequent section of this report for comparison with calcu- 
lation. 
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Table 1 

Operating Conditions NASA Langley Tests 



22 X 5.5 

Tires 



Load 

lb 

Pressure 

psi 

Speed 

MPH 


2700 

83 

20 


3270 

83 

20 

NASA 

4000 

83 

20 

8PR 

4000 

115 

20 


4000 

150 

20 


2560 

55 

20 


3150 

55 

20 

NASA 

4000 

55 

20 

12PR 

4000 

110 

20 


4000 

78 

20 


4000 

78 

35 


4000 

78 

50 


Tab le 2 

Operating Conditions - University of Michigan Tests 

22 X 5.5 Tires 



Load 

lb. 

Pressure 

psi 



4350 

125 

20 

.UM 

3263 

87 

20 

8PR 

4350 

156 

20 


4350 

94 

20 


7100 

209 

20 

UM 

5325 

132 

20 

12PR 

7100 

261 

20 
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Table 3 

Operating Conditions - University of Mich igan Tests 


40 X 14 Tires 
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CALCULATION VS EXPERIMENT 


Calculated temperature distributions in the tire were first 
compared with data measured from both the NASA tests and the 
University of Michigan tests. In order to carry out the compu- 
tations, several material characteristics had to be estimated 
from the literature. These characteristics are given in Table 
4. 

Using a time increment of ten seconds, the computation was 
carried out for the temperature rise in these tires under condi- 
tions experienced by both the NASA tests and the tests conducted 
by the University of Michigan. Comparison between calculation 
and experiment is shown in Figures 19 through 24 for the thermo- 
couple locations A through F , where the experimental data are 
plotted as a function of time. The computations are carried 
out using equation 39 and the material constants given in Table 
4, along with the specific geometry for these tires. 

It should be emphasized that while the outside dimensions of 
these two ply-rating tires are the same, their internal charac- 
teristics are quite different and hence the modeling is complete- 
ly different. The elements representing the heat transfer cells 
are not the same, nor is the nature and extent of the carcass 
area. The 8PR version of the tire has only one bead, while the 
12PR version has two beads. For these reasons the two sets of 
comparisons should be considered as being taken on two separate 
tires . 

Similar comparisons between calculation and theory are 
given for the two 40 x 14 tires, and again plots of temperature 
rise are used as a basis. These data are presented in figures 
25 through 36. The material constants used for the 40 x 14 
tire are also given in Table' '4 and once more some of the input 
data must be estimated from literature sources. 
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Table 4 

Material Characteristics of 
22 X 5.5 and 40 x 14 Aircraft Tires 


Region* 

■ 

■ 

wm 

IIIIQffiR 


K 

B 


Tread (II) 


335 

0, 15 

— 

-4 

5x10 ^ 

1.0 

0.5 

Carcass 

(I, III) 

(1) 


0.15 

(2) 

-4 

5x10 

1.0 

0.5 

Sidewal 1 
Rubber (IV) 


335 

0.15 



1.0 

0.5 

Bead (V) 

15x10® 


0.03 






* See fig. 6 

(1) See Appendix E. For shear modulus G of rubber use 335 psi. 

(2) Calculated by cosine law using a = 35° crown angle. 

(3) Values of tan 6^ are given for 25°C, For higher temperatures, 
tan6 is calculated by the approximate expression 


tan 6 = tan 6 

o 


-.O1A0 

e 


where A0 is the temperature rise in 


above 


ambient . 


E = Young’s modulus G = Shear modulus of a = Cord 

of rubber, psi rubber, psi angle, 

K = Thermal p = Density C = Specific 

conductivity , 2 ^ heat 

cal -cm/^C-cm'^-sec 


deg. 
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It should be emphasized that some of the input data needed 
for these computations is not readily available to the casual 
user of this program. In particular this holds true for the quan- 
tity tan 6 , which directly controls the rate of heat generation. 

This quantity is not well known for cord-rubber composites, and 
in this case it has been estimated from test data on similar com- 
posites. Both the strain and temperature dependence of tan 6 are 
also important. We have attempted to estimate the temperature 
dependence term from very sketchy published data. The strain 
dependence is left undefined, but known to be a large effect. 

These shortcomings point up the need for continued research in 
material properties of cord-rubber composites. 

Bending stiffness of the composite sidewall of the tire is 
important, particularly in the area near the bead. It is very 
difficult to calculate this quantity with any accuracy since it 
depends strongly on the location of the neutral axis, which, in 
turn, is influenced by the state of tension in the cord-structure. 
Experiments show bending stiff nesses which are manytimes smaller than 
those which would be predicted by simple beam or plate calcula- 
tions based on properties measured in tension. For this reason 
care should be taken in the bead region of the tire to insert 
proper bending stiffnesses. 

Runway friction influences to some extent the membrane 
stresses in the center of the tire contact patch. This is very 
difficult to estimate quantitatively. It is discussed analy- 
tically in Appendix J, but the actual friction levels available 
in the tests which are quoted can only be estimated in calcula- 
tion. 

As can be seen from figures 19 through 36, the calculations 
generally agree well with the experimental data over a wide 
range of inflation pressures and tire deflections. In parti- 
cular the calculations predict the ranking of hottest areas 
in the tire at a particular location, i.e., from outside to 
inside and vice-versa. Further, the calculations support the 
reversal of inside-outside temperature gradient between sidewall 
and tread, a phenomenon observed experimentally. 
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From a study of the comparisons between these computations 
and the experimental data^ it is concluded that the major trends 
in the temperature rise problem can be calculated adequately for 
relatively short term taxi-takeoff conditions. More emphasis 
must be placed on obtaining good thermal properties of the mat- 
erials and better heat transfer coefficients in order to improve 
the accuracy of the predictions. 
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Fig. 19, Measured and calculated tire temperature profiles 


TEMPERATURE RISE 


HRE:22x5.5-8PR O INSIDE 

SPEED:20MPH A MIDLINE 

FREE ROLLING □ OUTSIDE 

UM TEST CALCULATED 


SURFACE: I20IN. DIA.DRUM 
Fz : 4350/4350 = f00% RATED 
Sz : 1.341 / 1. 344 = 99.8 % rated 
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Fig. 20, Measured and calculated tire temperature profiles 








TEMPERATURE RISE 


TIRE:22x5.5-8PR 
SPEED: 20MPH 
FREE ROLLING 
UM TEST 


O INSIDE 
A MIDLINE 
□ OUTSIDE 
— CALCULATED 


SURFACE: 120 IN. DIA. DRUM 
Fz : 4350/4350 = 100% RATED 
Sz : 1.656/1.344= 123% RATED 
Po : 94/115 =82% RATED 
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Fig, 21. Measured and calculated tire temperature profiles 




TIRE 5 22x5.5 12 PR 
SPEED:20MPH 
FREE ROLLING 
UMTEST 


O INSIDE 
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□ OUTSIDE 
CALCULATED 


SURFACE : 120 IN STEEL DRUM 
Fz= 7100/7100 100% RATED 
Sz 304/1. 312 99%RATED 
Po* 209/190 110% RATED 
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Fig. 22. Measured and calculated tire temperature profiles 



TEMPERATURE RISE 


TIRE:22x5.5-!2PR O INSIDE SURFACE: 120 IN. STEEL DRUM 

SPEED: 20MPH ^ MIDLINE F, : 7100/7100 = 100% RATED 

FREE ROLLING □ OUTSIDE : 1.134/ 1.312 =86% RATED 

UM TEST — CALCULATED Pq : 261/ 190= 137% RATED 
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Fig. 23, Measured and calculated tire temperature profiles 
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TEMPERATURE RISE, 


TIRE = 22x5.5-12 PR ° INSIDE 

SPEED! 20 MPH ^ MIDLINE 

FREE ROLLING □ OUTSIDE 

UM TEST CALCULATED 


SURFACE ! 120 IN STEEL DRUM 
fz ■■ 5325/7100 =95% RATED 
Sz ' 1.306/1.312 = 99.5% RATED 
Po ! 132/190 =69.5% RATED 
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Fig. 24, Measured and calculated tire temperature profiles 
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TEMPERATURE RISE 


TIRE:40xl4-22PR O INSIDE 


SURFACE: 120 IN. DIA. DRUM 


SPEED: 20 MPH 
FREE ROLLING 
UM TEST 


A MIDLINE 
□ OUTSIDE 
— CALCULATED 


Fz ; 23,500 / 25,000 = 94 % RATED 
Sz : 3.378 / 3.378 = 100 % RATED 
Po :l60/l55 = l03 7« RATED 
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Fig. 25, Measured and calculated tire temperature profiles 





TEMPERATURE RISE 


TIRE:40xl4-22PR 
SPEED: 20 MPH 
FREE ROLLING 
UM TEST 


O INSIDE 
^ MIDLINE 
□ OUTSIDE 
— CALCULATED 


SURFACE: 120 IN. DIA. DRUM 
Fz ; 21,000 / 25,000 = 84 % RATED 
Sz : 3.378 / 3.378 = 100 % RATED 
Pq : 145/155 = 94% RATED 
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Fig. 26. Measured and calculated tire temperature profiles 


TEMPERATURE RISE 


TlRE:40xl4-22PR O INSIDE 


SURFACE: 120 IN. DIA.DRUM 


SPEED: 20 MPH 
FREE ROLLING 
UM TEST 


^ MIDLINE 
□ OUTSIDE 
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Fz ; 25,000 / 25,000 = 100 % RATED 
8z : 3.378 / 3.378 = 100 % RATED 
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Fig. 27. Measured and calculated tire temperature profiles 





TEMPERATURE RISE 


TIRE: 40x14-22 PR O INSIDE 
SPEED: 20 MPH ^ MIDLINE 
FREE ROLLING □ OUTSIDE 
UM TEST — CALCULATED 


SURFACE; 120 IN. DIA. DRUM 
Fz ; 29,300 / 25,000 =117% RATED 
Sz : 4.054/3.378 = 120 % RATED 
Pq *. 160 / 155 = 103% RATED 
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Fig. 28, Measured and calculated tire temperature profiles 










TEMPERATURE RISE, °F 


TIRE 40x14 -22PR ° INSIDE SURFACE: 120 IN. DIA, DRUM 

SPEED; 20 MPH A MIDLINE : 35,400/25000 = 142% RATED 

FREE ROLLING o OUTSIDE 8, ; 4.729/3.378= 140% RATED 

U M TEST — CALCULATED A . (60/155= 103% RATED 



Fig. 29. Measured and calculated tire temperature profiles 
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TEMPERATURE RISE, 


TIRE: 40x14-22 PR 
SPEED: 40 MPH 
FREE ROLLING 
UM TEST 


o inside 

^ MIDLINE 
a OUTSIDE 
— CALCULATED 


SURFACE: 120 IN. OIA. DRUM 
Fz • 25,000 / 25,000 = 100 % RATED 
Sz . 3.378 /-3.3T8 = 100 % RATED 
P(j : 175/ 155= 113 % RATED 
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Fig. 30. Measured and calculated tire temperature profiles 








TEMPERATURE RISE 


TIRE:40kI4-28PR O INSIDE SURFACE; 120 IN. DIA. DRUM 

SPEED: 20 MPH ^ MIDLINE Fz ; 30,800/33,500 = 92% RATED 

FREE ROLLING ° OUTSIDE Sz : 3.459/ 3.484 = 99 % RATED 

UM TEST —CALCULATED Pq : 207/200= 104 % RATED 
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Fig. 31. Measured and calculated tire temperature profiles 




TEMPERATURE RISE 


TIRE:40xl4-28PR 
SPEED: 20 MPH 
FREE ROLLING 
UM TEST 


O INSIDE 
^ MIDLINE 
□ OUTSIDE 
— CALCULATED 


SURFACE: 120 IN. DIA. DRUM 
Fz : 27,800/33,500 = 83% RATED 
8z : 3.484 / 3.484 = 100 % RATED 
pQ : 187/200 = 94 % RATED 
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TIME, SEC 


Fig, 32. Measured and calculated tire temperature profiles 
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TEMPERATURE RISE 


TIRE;40xl4-28PR O INSIDE 

SPEED: 20 MPH ^ MIDLINE 

FREE ROLLING □ OUTSIDE 

UM TEST — CALCULATED 


SURFACE; 120 IN. DIA. DRUM 
Fz ; 33,500/33^00= 100% RATED 
Sz : 3.484/3.484 = 100% RATED 
Pq ; 227/200= 114% RATED 
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Fig. 33. Measured and calculated tire temperature profiles 



TEMPERATURE RISE 



TIRE; 40xl4-28PR O INSIDE SURFACE; 120 IN. DIA. DRUM 

SPEED; 20 MPH A MIDLINE Fz 38,200/ 33,500 = 114 % RATED 

FREE ROLLING D OUTSIDE 8z ; 4.151/ 3.484 = 119 % RATED 

U M TEST — CALCULATED Pq •• 207/ 200 = 104 % RATED 














TEMPERATURE RISE 



TIRE: 40xl4-28PR O INSIDE SURFACE: 120 IN. DIA. DRUM 

SPEED: 20 MPH ^ MIDLINE Fz 46,300 / 33,500 = 138 7o RATED 

FREE ROLLING □ OUTSIDE 8z : 4.877/ 3.484 = 140 7o RATED 

U M TEST — CALCULATED P© : 207/ 200 = 104 7o RATED 
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Fig. 35. Measured and calculated tire temperature profiles 



TEMPERATURE RISEV 


TIRE; 40x14-28 PR O INSIDE a SURFACE; 120 IN. DIA. DRUM 
SPEED: 40 MPH ^ MIDLINE ° Fz ; 33,500/ 33,500 = 100% RATED 

FREE ROLLING Q □ OUTSIDE O 8z : 3.484 / 3.484 = 100 % RATED 

UM TEST — CALCULATED Po ; 227/ 200 = 114 % RATED ^ 



TIME, SEC 


Fig. 36. Measured and calculated tire temperature profiles 
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CONCLUSIONS 

A method has been developed for approximating the internal 
temperature distribution in an aircraft tire due to unyawed free 
rolling under load. The method is based on an approximate stress 
analysis of a material point in the tire as it rolls through the 
contact patch. From this analysis the mechanical work done on 
the material point is converted into the heat generated due to 
the loss characteristics of the materials at that point. Thus, 
each material point of the tire has its own cyclic heat generation 
each time it rolls through the contact patch and is stressed. 
Furthermore, due to the high frequency of cylic stressing compared 
with the relatively slow diffusion times of heat in polymeric 
materials, the temperature distribution in the tire is considered 
to be axisymmetric . 

With a rate of heat generation established at each material 
point on the cross section of the tire, heat diffusion through- 
out the carcass is calculated by the usual numerical methods which 
will involve loss of heat in the tire through various mechanisms 
such as convection at the outer and inner surfaces, conduction 
through the bead area into the flange of the wheel, and conduction 
through the tread surface into the runway material. All of these 
mechanisms modify the normal diffusion process so that at any time 
tire internal temperatures are dependent not only upon internal 
rates of heat generation but also on the appropriate loss charac- 
teristics at the tire surfaces. 
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An extensive experimental program was carried out jointly 
between the University of Michigan and the Impact Dynamics Branch, 
Langley Research Center, NASA. This program consisted of tempera- 
ture measurements made at internal points in 22 x 5,5 and 40 x 14 
size tires, each with two different ply ratings. Due to differences 
in dimensions, section thickness and in distribution of cords, 
these four tires each require a completely separate geometric and 
material property description. Some of the experimental data were 
taken at the NASA Wallops Flight Center using the NASA tire test 
truck, operating on asphaltic and concrete runway surfaces. The 
University of Michigan data were taken on the 120-inch steel dynamo- 
meter drum at Wright Patterson AFB, Ohio, where higher tire loads 
could be obtained. The data presented in this report are mostly 
taken from these latter tests. 

All data were obtained from thermocouples buried in the tire 
carcass prior to retreading the tires. These thermocouples were 
connected to conventional recording apparatus , and gave continuous 
readings of tire temperature at selected internal points. 

Experimental data from these measurements are compared with 
data calculated over a wide range of tire loads, inflation pres- 
sures, speeds, and deflections. Generally speaking, the analyti- 
cal model gives good agreement with measured values when the tire 
is operating in the range of its normal rated load, pressure, and 


deflection. 
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APPENDIX A 

Energy Loss as Related to Strain Energy 

The question of the formulation of the hysteretic loss in 
materials is not well defined in the literature. Elementary 
considerations, such as discussed in the body of this report, 
imply that hysteretic loss may be expressed as a fraction of 
strain energy. In the case where the hysteretic loss is caused 
by cycling between two non-zero stress states, it is not immed- 
iately clear whether the appropriate strain energy is the 0-1-2-4 area 
or the area 1-2-3 shown in Fig. A1 . 



Fig. Al. Idealized stress strain curve 

In order to resolve this question, loss characteristics of 
a typical bias cord-rubber composite were measured over a range 
of prestrains using a constant cyclic strain excursion. If the 
strain energy upon which loss is based is the area 0-1-2-4 of 
Fig. Al, then the hysteretic loss should increase with prestrain. 
On the other hand, if the area 1-2-3 defines the appropriate 
strain energy then the hysteretic loss should be independent of 
prestrain. The results of the experimental program showed a 
constant tan6 over a range of prestrains, thus validating the 
concept that the appropriate area is 1-2-3. Details of the ex- 
periment are given in Table Ai, 
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Table A1 

Details of strain energy experiments 
Test A B C D E 


Prestrain (in) 

0.2 

0,4 

0.6 

0.8 

1.0 

Stroke (in) 

+0.1 

+0. 1 

+0.1 

+0.1 

+0.1 

tan 6 

0.151 

0.147 

0.148 

0. 144 

0.139 


Specimen: 1000/2 Polyester cord in Rubber Cylinder 2.6 in. 

diameter x 15 in. long. Frequency 10 Hz. 

Cord angle 45^^. Loading axial. 
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APPENDIX B 


Cross-Section Curvatures 

As discussed in the body of the report, the inflated shape 
of the tire cross-section must be obtained for the purpose of 
determining the geometry of the midline of the tire carcass. 

This geometry enters into the computation of the membrane 
stresses . 

Several different methods were examined for determining the 
cross-section shape. Since aircraft tires in general are of bias 
construction, the Hofferberth equations of reference 4 were exa-^ 
mined for their accuracy of prediction. It was found that they 
did an adequate job of predicting the inflated shape of the tire 
in its upper half, in the vicinity of the tread and shoulder region, 
but seem to fall quite short of good accuracy near the lower side- 
wall and bead area where rim and bead constraint effects were 
greater. Deviations in these regions were large enough so that 
the Hofferberth equations were deemed inadequate. 

An examination of the differences in such tire dimensions as 
overall diameter and section width between the uninflated and 
inflated tires showed that these were sufficiently small to the 
extent that to a good approximation the tire cross-sectional shape 
when inflated, could be approximated from its uninflated cross- 
section. Using sections cut from the tire and making appropriate 
changes for the effect of inflation, the radial positions of the 
various element centroids could be approximated, as could the 
tangent angles needed for computation such as shown in Fig. 3. 

Such drawings were made for both types of 22 x 5.5 and 40 x 
14 tires. Examination of these drawings showed that the sidewall 
region can be approximated rather well by a single circular arc 
beginning at the bead and ending at the edge of the tread. The 
tire carcass centerline in the tread region can also be fit by 
another circular arc segment blending into that used for the side- 
wall. Such a representation has some advantages since it allows 
rapid calculation of membrane stress levels. 
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Calculation of the tire deformed shape requires considerably 
more effort. In the case of the tread region, the tire cross- 
section is considered to be flat and in complete contact with the 
runway surface. However, the sidewall region is sever ly deformed. 

The point of demarcation between the tread and sidewall region 
is clearly defined in most aircraft tires by the molded pattern 
of the tread elements. In the cases of all four tires used for 
this study this point is clearly evident, allowing it to be 
marked . 

The region outside of contact deforms substantially, and 
one method of obtaining an accurate picture of the geometry of 
this deformation is the use of two-dimensional finite element 
programs suitable for handling large deformation processes. 

Again one of the shortcomings in such a use is lack of good 
elastic property data on the composite carcass stiffness. 

A simpler computation scheme can also be employed. The 
tire sidewall can be considered to be a curved arch as shown in 
Fig. Bl, and loads can be applied to the upper end of the arch 
sufficient to cause the known deflections which must take place 
there, as determined by the deflection imposed on the tire. 

Fig. Bl illustrates the geometric variables involved here. 



Fig. Bl. Geometry of curved arch 
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If the bending stiffness of the arch is known at the various 
points throughout its length, then linear energy analysis should 
be sufficient to determine the values of the loads needed for 
this deflection. From these loads, changes in radii of curva- 
ture can be obtained, as well as deflections of the intermediate 
points. By a similar analysis, the tangent angle of the deformed 
shape can also be obtained. Such an analysis was carried out for 
both the 22 X 5.5 and the 40 X 14 tires, and these shapes showed 
relatively good agreement with the known section widths of these 
tires at the center line of their deflected contact patch. 

In both cases it was necessary to modify the usual calculation 
somewhat. The linear computation of such large deflections re- 
sults in an unexceptably large change in arc length of the curved 
beam. For that reason, the vertical deflections were calculated 
by the usual linear method, and the horizontal deflections were 
obtained by assuming that the arc lengths between each of the 
elemental points were held constant. A typical cross-section 
drawing for a 40 X 14 tire in its undeformed and deformed posi- 
tion is shown in Fig. B2. 

Finally the simplest method for representing the sidewall 
shape is to consider the relatively thin sidewall as a membrane 
enclosing a two dimensional opening extending from the bead to 
the edge of the tread region in contact with the runway. This 
membrane will form a circular arc whose radius can be determined 
from the length using an average curvature, which can then be 
used to estimate both bending and membrane effects as tire de- 
flection or pressure vary. 

If one uses the circular arc model for both deformed and unde- 
formed tire shapes, then much of the needed tire geometry can be 
calculated from the arc properties. 

It is assumed that the geometry of the tire carcass can be 
described by two intersecting circular arcs. For the undeformed 
tire these arcs describe the carcass centerline from the crown 
of the tire to the shoulder, R , and from the shoulder to the 
bead, R^, as shown in Figure B3, 

It is assumed that the x-y coordinates of the shoulder 

bead (x , y ) , the y-coordinate of the carcass at the 

O O O £5 

crown, y^, and the length of the carcass centerline from the 
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shoulder to the bead, S^, are known. All of these geometric 
properties can be measured from appropriate tire cross-section 
profiles. The location of the shoulder is somewhat arbitrary. 
However, for this analysis was taken to be 0.75 of the section 
half-width w, and y to be the intersection of x=x and the car- 
cass centerline. 



Fig. B3, Tire cross-section illustrating 

idealized curvatures R and R- 

c 1 
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From geometry: 


2 j ^ 2 

X + (y -y ) 
c -‘^m 

2(y -y ) 

■^m -^c 


(Bl) 


. -If c. 
sin i:^) 


(B2) 


S = R *6 
c c c 


(B3) 


To find the radius of curvature of the undeflected sidewall 
one must first determine H; 


H = /(x^-Xg)|(y_,-yg)2 


(B4) 


Noting that = Sj^/0 and that H/2 - R^ sin (-|-) , the angle 0 can 
be found from: 


ft 

0 sin(J) = 0 


(B5) 


Solving (B5) for 0 then allows the calculation of R^, 

To calculate the radius of the deflected sidewall it is 

first assumed that the deflection is greater than y -y and that 

the carcass changes from a circular arc of R^ to a flat surface 

of half-length S . It is thus assumed that the x-coordinate of 

c 

the shoulder in the deflected position is and the y-coordinate 

is (y -A) . With the new coordinates of x and y a new H can be 
■^m c 

calculated from (B4). Also from the figure: 


^f ^f 

R^ sin(^) = 


(B6) 
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Substituting from (B7) into (B6) : 


(B7) 


sin 




^B8) 


The solution of this transcendental equation for gives the 
radius of curvature of the deflected sidewall. 
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APPENDIX C 


Bead Tension Analysis 

Consider the sidewall of an aircraft tire as shown in fig. 
Cl, first in the undeformed and then in the deformed state. 





Fig. Cl. Sidewall geometry used in bead tension analysis 

The membrane forces in the (p direction are given by Eq. (2) 
as 


= 

B 


„ 2,2 

P ,z, -b. V 

o ( b 1 ) 

2 z, sin3. 
h 1 


= 

^ B 


r, 2,2 

P ,z, -b^ , 

o ( b f ) 

2 Zj^sing^ 


The angles 3 are obtained from the tire geometry analysis 
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Bead tension may be split into radial and axial components 
as in fig. C2. The axial component is carried by the wheel 
flange. The radial component causes bead tension in the amount 

= (N^ sin3) (C4) 

using a uniformly loaded ring as a model of the bead. 



Fig, C2 . Details of bead area 

Denote the gross area of the bead wires as and the modu- 

lus of the bead bundle based on this gross area as E . Then the 

B 

strain energy stored during one half cycle is 

f 1 
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APPENDIX D 


Sidewall Membrane Shear 

The membrane forces and vary as the tire rolls through 

the contact patch. The variation in Ng causes a shear stress 

to be required for equilibrium purposes. This shear stress gives 

rise to an energy storage, hysteretic loss and heat generation 

process which is similar to those found by other mechanisms. 

Consider the tire cut at its centerline and spread out 

flat as shown in fig. Dl. Equilibrium of forces requires the 

existence of a membrane shear force N as shown. This is caused 

xy 

by the variation in N with circumferential position. 


/(Ng+ANg)ds 


Nxydx 


dx 

Rdd 


ds 


— ►/N^ds 

-t 


Fig. Dl. Idealized view of tire section 


Next consider the relation between Ng and N^ . For a net 
structure in which loads are carried by the cords alone, Fig. D2 
illustrates the division of cord load between N^ and Ng . If T 
represents a typical cord load, it may be shown that 



(Dl) 


Fig. D2 


Idealized net structure 




81 


Since is calculated as part of the contribution to strain 
energy in the sidewall regions I and III, then it is available 
for computation of . At any location the variation of Nq 
follows the variation of , which is approximately of the form, 
shown in fig. D3. 



b + 


Fig. D3. 


Variation of N. with circumferential position 
6 


Letting the variation of with respect to circumferential posi- 
tion X be approximated by 


2 2 

= a + bd-e"*^ ^ ) 

u 


(D2) 


where 




and c = constant to be determined later. 

Since the maximum value of N depends on the maximum gra- 

xy 

dient of N^, then it can be determined by differentiation of 
Equation D2. 


dN^ 

dST 


„ 2 2 
^,2 -c X 
2bc xe 


d^N 


dx 


0 ov.2n..22 22 

= 2bc [l-2c X , “C X 
2 ] e 


The maximum of dN„ occurs at x - 


9 


dx 


c/7. 
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and results in 


dN 

( ~) 

dx max 


2b - — e“^ = 0.858 be 

c/2 

0.858c[ (NQ)^-(NQ)g] cot^ a 


From plane equilibrium. 


(D3) 


dN„ dN 

6 _ xy 

dx ~ dy 
so that 


(D4) 


dN 


xy 


0.858c[(N^)j^-(N^)g] cof 


a*dy, where y 


is measured along the circumference of the tire meridian. Cot a 
varies with this position, so that the net value of N^^ at any 
computational cell is given by 

N 


xy 


= 0.858c[(N^)^-(N^)3] 


r 


cot a • dy 


(D5) 


We now approximate the constant c representing the charac- 
teristic distance needed for change in the tensile forces in the 
circumferential direction, assuming that 


c — 


(one contact patch length) 


-1 


Since the contact patch length is a function of the tire 
size and tire deflection, a separate analysis of it is carried 
out in Appendix F, where it is shown that 

i 

contact patch length — [3A(D-A)]^ 

v^here A = tire deflection 
D = tire diameter 
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From this the strain energy can be formed 


U 


1 /h^G 

2 xy' xy 


(D6) 
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APPENDIX E 


Aircraft tire construction is entirely bias ply and the ply 
angle range is limited between 30^ and 60° for half angles ct as 
illustrated in fig. 7. Previous studies by the authors Iref. 2 ] 
and by Walter [ 3] have shown that within this range of angles 
the tensile stiffness of the composite bias ply structure can be 
represented by a simplified analytical model using cords which 
are inextensible but which are encased in an elastic matrix. 

This is often given the name "net theory" . 

One of the shortcomings of such a theory is that some of 
the other elastic constants needed for this analysis are not 
predicted well by use of net theory. In an effort to improve 
the accuracy of those computations a somewhat more advanced 
theory was constructed around the general framework of a net 
analysis. In this theory, however, the cord is now assumed to 
have a certain amount of extensibility so that in extreme cases 
the stiffness characteristics will not approach indefinitely 
large values such as is possible with a pure net theory. 

Coordinate directions and cord angles are given in fig. El. 



Fig, El. 


Detail of idealized bias 
net structure 
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Let a = cord half angle with respect to 0, and let 
V E (1-V ) 

4) = ^ = stiffening Parameter 


where 


V = volume fraction of cord 

c 

E = cord modulus 

c 

G = rubber modulus 


Then 


X = 


(tan^g+1) (2tan^g-l) 
2 (tan^a+1) ^ + (|) 


and the composite (El) 


elastic modulus in the 4 direction is found from 


2G 


- [ 


I (1+tan^a)^ + 2x(l-tan^a) 


+ (1+tan^a) ^ 


+ I (l-X)^(l+tan^a)^ 


+ X^(|) 


(E2 ) 


The Poisson's ratio is 

Pq, = (1+tan^a) (1-X)”1 

0 cp 


(E3) 


by 


The angle distribution in a bias ply aircraft tire is given 


cos a = — cos a 


where 
a _ = 


crown angle 
crown radius 


a 


Fig. E2. Tire cross section 
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APPENDIX F 

Referring to fig, Fl/ the length of the contact may be app- 
roximated as follows 



R (1-cos 6 ) = A 

cos 0 = 1 - 4- 


Sin^0 


2 ^ ^ 
^ ' r2 


1 

2 


Rsin0 



4A(D-A) 


Frcm reference 5 , the empirical expression for the con- 
tact patch length L is given as o*85 times the geometric inter- 
section length t, so that 


= (0.85)^x4 (D-A) =3A(D-A) 


from which 


L 


/3A(D-A) 
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APPENDIX G 


Independent of other effects, the tire cycles between zero 
tread compressive stress and a maximum compressive stress. 

Measurements indicate a maximum pressure/ unit length at the 
centerline of contact, approximating 1.25 x average inflation 
pressure x width. 

Examine three possible pressure distributions, all giving 
rise to the same load/unit length. See Fig. G^-l. 


Stored energy 

K 


w/2 

u 


p dy assuming E = constant across the 

tread . 


(a) E 


I 


^ 2 _ 2 

^ j (1.25p^) dy - 1.5625p^w for constant pressure 


(Gl) 


(b) Parabola 1; p = p [1“(-^)^] 


m w 


r w/2 

J ’ 

•'-w/2 


p dy = pm iy 


2f 16y^ | 2 . 8 

' = pm ) 

3w 5w J 

-w/2 


But 


/ 


w/2 


pdy =1.25 p^w 


-w/2 


rw/2 

-w/2 


w/2 . 2 

p [1 - dy 

w 


2p w 


Pm = 7 ^ P. 


(G3) 


2 8 

p X -s-?W 
TT, 15 

m 


9 25 2 8 

4 16 ^o 15 


= 1.875p w = 1.2 X E 

^o a 


(G4) 
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(c) Parabola 2: P=Pj^(^^) 

w/2 
-w/2 


j ' 

•'-w/2 


pdy = 1.25 


p w 


= P. 


m 


t 


A 2 
2 


dy = 


m 


3w^ 


E 


c 





2.8125 


2 

p w 


1.80 E 

a 


We adopt a conservative value of ^ ‘ ^^ 2 ^ * “ 1 « 5 x uniform 

case. Assume uniform contact pressure of 1.25 x p^. Then 

multiply by 1.5 factor for shape distribution, giving a total 

value of 1.25x1.50 p = 1.875 p . 

^o o 



<L 



Fig, Gl. Transverse pressure distributions 
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APPENDIX H 


Thermal Analysis 


In the thermal analysis , account is taken of 
conduction to neighboring elements, conduction to the road 
surface, convection to the outside air and to the inside air, 
and heat generation within each element. In writing the heat 
balance equation for each element it is necessary to provide 
effective heat transfer and heat generation coefficients for each 
element. The necessary coefficients are: 


AV, AH 


a- At 

Cjj-Pjj-vd 




93 


where : 

AV, AH 

QB 

HF 

HR 

HA 

HC 


V 

a 


d 

At 


4 


K 

K 

C 

C 

p 

p 

h 


R 

A 

R 

A 

R 

A 

R 


h 


A 


h 


c 


h 


f 


HF 


Kjj.a 


K . . a 
A 




C^PA'Vd 


hj-a 

Spr-'^ 


- effective heat 
elements 

- effective heat 
effective heat 
interface 

- effective heat 
face 

- effective heat 
air interface 
effective heat 
air interface 


conduction coef f iceints , neighboring 

generation coefficient, each element 
transfer coefficient, tire-rim 

transfer coefficient, tire-road inter- 
transfer coefficient, tire-outside 
transfer coefficient, tire-contained 


= volume of each element 
= area across which heat is transferred 
= distance between nodes 

= time increment - selected 

= rate of heat generation each element - see earlier sections. 

2 

= thermal conductivity of rubber = 5x10""^ cal-cm/°C*cm -sec 

2 

= thermal conductivity of rim = .48 cal • cm/ "C* cm *sec 
- specific heat of rubber = 0.5 cal/gm °C 

= specific heat of rim = .23 cal/gm °C 

3 

= density of rubber = 1 gm/cm 

3 

= density of rim = 2.7 gm/cm 

= film heat transfer coefficient, tire-road interface 
= .25 cal/® C* cm^ • sec 

= film heat transfer coefficient, tire-outside air inter- 
face = (1.3x10 "^+4x10 *^V( cm/sec)) cal/ ®C • cm^ • sec 

= film heat transfer coefficient, tire-contained air 

-4 2 

interface = 1x10 cal/°C*cm *sec 
= film heat transfer coefficient, tire-rim interface 
= .5 cal/“C cm^*sec 
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Typical heat balance equations for the individual elements of 
the 22x5.5 8PR tire are listed below where each equation has been 
solved for the change in temperature which takes place during a 
given increment of time. 

The inside elements all have an equation of the form: 

ATj ^ = QB(I,1) + AH(I.l)(Tj 2 ~Tj + AV( 1 , 1 _ j^-Tj 3 ^) 

+ AV(I-1,1)(Tj_ 3^ ,-Tj ,) + HC(I)(TC-Tj. 3 ^) 

The middle elements all have an equation of the form: 


ATj. 2 = QB(I, 2 ) + AH(I, 2 )(Tj g-Tj. 2> + AH( 1 , 1 ) (Tj _ 

+ AV(I.2KT3^3_2-^I,2^ + AV( I-l , 2 ) g'Tj , 2 > 

The outside elements not in contact with the road surface or rim 
have an equation of the form: 

ATj 3 = QB(I, 3 ) + AH(I, 2 )(Tj 2 -Tj 3) + AV( I , 3 ) ( T^^^ 3~Tj 3 ) 

+ AV(I-l, 3 )(Tj_^ 3 ~Tj 3) + HA(I)(TA-Tj 3) 

The outside elements in contact with the rim have the same form 
except HA is replaced by HF and TA by TF. 

The outside elements in the tread region are assumed to be in 
contact with the air 90% of the time and in contact with the 
road surface 10% of the time. The heat balance equation for 
these elements is of the form: 


AT 


1,3 


QB(I, 3 ) + AH(I, 2 )(Tj 2 -Tj 3) + ^ ^ ^ ( Tj^^ 3-Tj 3 ) 

+ AV(I-1,3)(Tj_^^3-Tj 3 ) + 0.9 HA(I)(TA-Tj 3 ) 

+ 0.1 HR(I)(TR-Tj 3 ) 
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In this analysis the outside air temperature TA, the rim tempera- 
ture TF, and the road surface temperature TR are assumed to 
remain constant. However, the change in the contained air tem- 
perature TC is approximated as follows: 

h *A* At- (T . -TC) 

G T 

TC = 

c c 

where : 

= ,24 cal/gm° c 

3 

= .0013 gm/cm 

= weighted average of inside surface temperatures 

For the 22x5.5 8PR tire analyzed, the approximate values of A 

2 3 

and V are 30.4 cm and 90.35 cm , respectively. Thus, for this 
tire; 

ATC = . 1078 (T^-TC) • At 
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APPENDIX I 


Contained Air Temperature 

Consider the tire cross-section as shown in fig. II. Let 
6 be the area weighted average of the temperatures on the inner 
surface of the tire over the region 0-1-2. Let y be the con- 
tained air temperature. Then the equation for temperature rise 
of the contained air is 


C ‘P'V- 


dt 


hA ( 0-y) 


1 



Fig. II, Tire cross-section 


where 


h = 
A = 


c 

P 

P 

V 


average heat transfer coefficient 
surface area of convection 0-1-2 
specific heat of air 
air density 

volume of contained air 


Let hA/CppV 


3. Then equation II becomes 


= 66 


(ID 


( 12 ) 


In the initial stages of tire heating the tire inner surface 
temperatures rise almost linearly, so that 0 = c x t_ where c is 
a constant. In this case equation (12) has a solution 
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equation (I2) has a solution 

y = e - (13) 

This solution lags behind the tire inner surface temperature. 

For longer times when 0 is determined numerically, the 
contained air temperature y must be calculated by the methods 
previously described using a finite difference form of equation 
(II) . 
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APPENDIX J 


Tread Friction 


The tread 
mated in a two 
action of both 
illustrated in 


region in contact with the runway can be approxi- 
dimensional sense as an elastic member under the 
tensile forces and surface friction. This is 
figure Jl. 


To ^ H H ^ 

n7¥mTh' ’m^ 


T, 


////// 


-H 


Let u = displacement in x direction 

Fig. Jl, Idealized section 

Let the tension T = AE of tire tread with 

runway 

where A = cross-sectional area 
E — modulus of elasticity 

Let K = shear stiffness of the tread rubber, which is assumed not 
to slip on the runway. 

The equation governing deformation is 


d^u K » 

- u = 0 


dx 


2 AE 


whose solution is 


^ Xx . _ -Xx 
u = C^e + C^e 


where 


^ AE 


Using the boundary condition 


(Jl) 


(J2) 


x= 


1 

2 


T 

o 


I 
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One obtains 


T(x) 



{ e^^+e“^) 


tr 3) 


This expression may be used to estimate the tension in the 
tread across the width of the contact patch. The value of K 
can be estimated from the known thickness, modulus and shear 
stiffness of the tread. In general the tread tension is reduced 
across the width of contact. 



2. Government Accession No. 


1. Report No. 

NASA CR-3629 


4. Title and Subtitle 


HEAT GENERATION IN AIRCRAFT TIRES 
UNDER FREE ROLLING CONDITIONS 


7. Author(s) 

Samuel K, Clark and Richard N. Dodge 


9. Performing Organization Name and Address 

University of Michigan 
Department of Mechanical Engineering 
Ann Arbor, MI 48109 


12, Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


15. Supplementary Notes 

Langley technical monitor: John L. McCarty 


3. Recipient's Catalog No. 

5. Report Date 

October 1982 


6. Performing Organization Code 


8. Performing Organization Report No. 


10. Work Unit No. 


11. Contract or Grant No. 

NSG-1607 


13. Type of Report and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 


16. Abstract 


A method has been developed for calculating the internal temperature distribution 
in an aircraft tire while free rolling under load. The method uses an approxi- 
mate stress analysis of each point in the tire as it rolls through the contact 
patch, and from this stress change the mechanical work done on each volume 
element may be obtained and converted into a heat release rate through a knowledge 
of material characteristics. The tire cross-section is then considered as a body 
with internal heat generation, and the diffusion equation is solved numerically 
with appropriate boundary conditions of the wheel and runway surface. Comparison 
with data obtained with buried thermocouples in tires shows good agreement. 


17. Key Words (Suggested by Authorls)) 

Aircraft tires 

Tire heating 

Tire thermal analysis 


18. Distribution Statement 


Unclassified - Unlimited 


Subject Category 03 


19. Security Classif. (of this report) j 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price* 

Unclassified | 

Unclassified 

98 

A05 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA-Langley, 1982 

















